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Abstract 

We propose a class of equity models whose volatility, Levy measure, and killing rate all have local 
stochastic state-dependence. In this framework we find a closed form solution for the price of any 
European-style option. Additionally, for a certain sub-class of models, we find an exact expression for 
the induced implied volatility smile. To illustrate our framework, we perform specific computations for 
hybrid CEV/Levy model (which we christen as the "CEV" model). 
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1 Introduction 



A local volatility model is a stochastic volatility model in "which the volatility at of an asset X is a function of 
the present level of X. That is, at = a(Xt). Among local sto chastic vola tility models, perhaps the most well- 



known is the constant elasticity of variance (CEV) model of lCoxl (| 19T5T ) . An extension of the CEV mo del to 
dcfaultable assets (the Jump-to-Default CEV or JDCEV model) is derived in lCarr and Linetskvl ([20061 ) . One 
advantage of these models is that they admit closed-form pricing formulas (written as infinite series of special 
functions) for European options. However, other than a possible jump to default, neither of these models (or 
any other local stochastic volatility model) allows the under lying asset to experience jumps (though jum ps 



Mendoza-Arriaga. Carr. and Linetskvl ([20101 )) 



can be inherited through Levy subordination, as described in 

Exponential Levy models, on the other hand, do allow the underlying price process to experience jumps. 
Furthermore, in an expo nential Levy s et ting, Europea n option-prices can be computed quickly as generalized 



Fourier transforms (see 



Lewis! ( 20011) 



Liptonl ((2(302)) . However, exponential Levy models are spatially 



homogeneous; neither the volatility nor the jump-intensity not have any local dependence. Thus, exponential 
Levy models do not exhibit volatility clustering or the leverage effect. 
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In this paper, we introduce a class of local Levy models. Such models allow for the volatility, Levy 
measure, and killing intensity of an asset to have stochastic local dependence. Within this framework, we 
find a closed form solution for the price of any European style option and, for certain sub-class of models, 
we find a closed form solution for the induced implied volatility smile. 

The rest of this paper proceeds as follows: in section[2]we present our class of models and our assumptions. 
In section [3] we derive a formula for the price of any European option. In section U we provide a formula for 
the implied volatility smile induced by a certain sub-class of models within our framework. And in section [5] 
we perform specific computations for a CEV-like local Levy model. Numerical examples are provided at the 
conclusion of the text. An appendix with some mathematical background is also provided. Some concluding 
remarks are found in section [51 



2 Model and assumptions 

We assume a frictionless market, no arbitrage and take an equivalent martingale measure P chosen by the 
market on a complete filtered probability space (Q, 3, {3^, t > 0}, P). The filtration {3^, t > 0} represents the 
history of the market. All processes defined below live on this space. For simplicity we assume zero interest 
rates and no dividends. Thus, all traded assets are martingales. We consider a risky asset X = exp(y) with 
lifetime £ whose dynamics are given by 

1/2 



dY t =a{Y t )dt+(a 2 +ea 2 1 ri(Y t )} dW t + zdN t {Y t ^,dz), Y =y, 

a(Y t ) = (c + £c l7? (r t )) - ^(a 2 + e al^Yt)^ - j^(F t _ , cb) (e 2 - 1 - z) , 

C = inf|t>0:^ (c +sc 1 r){Y s j)ds = e\. 

Here, (ao, o,i,cq, c\, e) are non- negative constants, r\ > is a function in Cq°(R) (the set of smooth functions 
on R with compact support), W is a Brownian motion, £ is an independent exponentially distributed random 
variable with parameter one, and dN t (Y t ^ 1 dz) is a compensated Poisson random measure 

dN t {Y t -,dz) = dN t (Y t ^,dz) - v(Y t -,dz), u(Y t _,dz) = v (dz) + ei]{Xt-)vi{dz). 

We assume that the locally-dependent Levy measure u(y, dz) satisfies 



min(l,z )v(y,dz) < oo, / e z v(y,dz) < oo, / \z\v{y,dz) < oo, 

J\z\>l J\z\>l 

for all y € M. The first integr ability condition must be satisfied by all Levy measures. The second integr ability 
condition is needed to ensure E Xt < oo for all T £ R + . The last integrability condition can be removed. 
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However, keeping this condition vastly simplifies the formulas that follow, as it will allow us to replace the 
indicator function that usually appears in Levy integrals with the constant one (i.e. I/i 2 i<i} — > 1). We list 
the main features of our class of models below. 

1/2 

• The process Y experiences local stochastic volatility with volatility function a(y) — (% + e af^O/)) 

• Jumps in Y of size dz arrive as a Poisson process with an instantaneous state-dependent intensity of 
v(Y t -, dz). The state-dependent Levy measure has the decomposition v{y, dz) = v^idz) +eT](y)ui(dz), 
where i>q and v\ are Levy measures. Note that both the jump intensity and jump distribution can 
change depending on the value of y. 

• The process Y is killed with a state-dependent intensity of k(Y t ) = (co + e cirj(Y t )). When Y is killed 
it is sent to a cemetery state A, from which it cannot return. 

• The drift function a(y) is fixed by the volatility and killing rate so that X — exp(F) is a martingale. 

• Both negative infinity and infinity are natural boundarie s for the process Y accordin g to Feller's 



boundary classification for one-dimensional diffusions (see 
Thus, the state-space of Y is R U {A}. 



Borodin and Salminenl ([20021 ) pp. 14-15). 



3 Option Pricing 

We wish to find the time-zero value u e (t, y) of a European-style option with payoff h(Y t )I^ >t -j at time t > 0. 
Using risk-neutral pricing we have 

u e (t,y)=E y h(Y t )I K>t} , 

where the notation Ej, indicates expectation starting from y. The function u £ {t, y) satisfies the Kolmogorov 
backward equation 

(~d t +A £ )u £ = 0, 7/(0,2/) = %), (1) 

where A £ is the generator of the process Y. The domain of A £ is defined as the set of / for which the limit 
lim t _>.o i (E y f(Yt)Ii£ >t \ — /(y)I{^>o}) exists in the strong sense. For any / £ C^°(R) the generator A e has 
the explicit representation 

A z = A + erjA u 

A i = \a i i (d 2 -d)+Ci{d-l)- [ Ui(dz)(e z -1- z)d+ [ v^dz) (e zd - 1 - zd) , *G{1,2}, 
domGAi) = C °°(R). 
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The symbol d (without a subscript t) indicates a derivative taken with respect to y. Note that the operator 
d is the generator of the translations: e zd f(y) — f(y + z). 

Remark 1. Let "K be the Hilbert space L 2 (M.,dy) with inner product (f,g) := L f(y) g(y) dy and norm 
11/11 = v ifi /)• Note that the operators {Ai,i — 1,2} are normal operators in 3~C and satisfy the following 
(improper) eigenvalue equations (neither Ai have any proper eigenvalues) 

A ip\ = (f>\ip\, Axif;\ = x\i>\, 

where 

4>\ = \al (-A 2 - iX) + co(iX - 1) - / ^o(^)(e 2 - 1 
1 Jr v 

Xa = £a? (- A 2 - iX) + ci (*A - 1) - / (dz) (e z - 1 

The eigenfunctions satisfy (-0A, Vv) = <5(A — /i). Note that Borel-measurable functions of normal operators 
(e.g., g(Ao)) are well-defined, as explained in Appendix |A"1 

We seek a solution to Cauchy problem (p} of the form 

oo 

U e = J2^ U n- (2) 

n=0 

Justification for this expansion will be given in Theorem [5] Inserting the expansion (J2|) into Cauchy problem 
(HJ) and collecting terms of like powers of e we find 

0(1) : (-d t + A )u = 0, u (0, y) = h(y), 

0(e n ) : {-d t + A )u n = -rjAm^, u n {0, y) = 0. 

The solution to the above equations is 

0(1): u {t,y) = e tA °h{y), 

0(e n ): u n {t,y)= [ dse^-^rt^Aiu^&y). 

Jo 

Next, using the equation (jlll) from appendix lAl we obtain 

0(1): u (t,y)= f dXe^i^h^xiv), 

0(e n ): Un (t,y)= [ [ dsd^e^ 8 ^ (^,r,Axu n ^(s,-))^(y), 

Jo Js, 



4>> 



/2tt 



z^jiX + / ua(dz) (^e 



z^jiX + / v\{dz)(^t 



1 — iXz 



lXz -1-iXi 
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After a bit of algebra, we find an explicit representation for u n (t,y) 

Un( *' y)= /W(n/ Afe ) (SrcX (n^^,^^))^^)^. (3) 

n+l 

We have now obtained a formal expansion ([I])-© for the price of a European option. The following theorem 
provides conditions under which the expansion is guaranteed to be valid. 

Theorem 2 (Option Price). Suppose e satisfies 

e2 < M W | 

- A€» ||,f . | & p 

/or some a > and 6 < 1. TTien the option price u E (t,y) is given by infinite series ([2]) where the {u n } are 
given by ((3|). 

Proof. See Appendix 151 □ 

Remark 3. Not that for certain <p\ and xa the bound (j4]) may be infinite (by choosing a large enough), in 
which case expansion ([2]) is valid for all e < oo. 



4 Implied Volatility 

Assumption 4. In this section only we assume v q = and Co = 0. Thus 

= \al (-A 2 - iA) . 

Note that, with this assumption, if e — 0, then X — exp(K) is a geometric Brownian motion. This will be 
key for our implied volatility analysis. 

We fix (t, y) and a call option payoff h{y) — (e v — e k ) + . Note that 

k— ik\ 

Wx,h)= ~ 6 Im(A)<-l. 

V27T \lX + \ z ) 

The following definitions will be useful: 
Definition 5. The Black- Scholes Price 

u bs . K + ^ R + is 

defined as 



,BS 



(a) := / dXe^W^h^x, <bf s (a) = l -o 2 {-\ 2 - 



Definition 6. The Implied Volatility is defined implicitly as the unique number o~ e S M + such that 

u BS {a e ) = u E , (5) 

where u e is as given in Theorem [5J 



Remark 7. Note that uo = u BS (ao). As shown in lLorid (|2012l) . when u e can be expanded as a power series 
whose first term corresponds to u BS , one can obtain the exact implied volatility surface corresponding to 



Remark 8. For < t < oo the existence and uniqueness of the implied volatility a e can be deduced by 
using the general arbitrage bounds for call prices and the monotonicity of u B . 

Remark 9. Note that u BS is an invertible analytic function that satisfies d p u BS (a) > for all a > 0. By 
the Lagrange inversion theorem, the inverse \u B ] _1 of such a function is also analytic. 

Clearly, u £ is an analytic function of e (we derived its po wer series expansion). It is a useful fact that the 



composition of two analytic functions is also analytic (see iBrown and Churchilll (jl996T ). section 24, p. 74). 
Thus, in light of Remark [9j we deduce that a e = [u ss ] _1 (u 6 ) is an analytic function and therefore has a 
power series expansion in e. We write this expansion as follows 



Taylor expanding u BS about the point Co we have 
u BS (a £ ) = u BS (a Q + 8 E ) 



5 £ = 5>V 

fe=l 



z — ' TV. 

oo 1 / oo \ ' 

71=1 ' \fc=l / 



71=0 

..BS 



d> BS (a ) 



,BS 



.BS 



BS 



oo -. 

71=1 

oo 

k=l 
oo 

M + E e * 



e e n°* 

oo / n 

eM e n°* 



ft> BS M 

u BS (a ) 



71=1 



fc=l 



oo 

n! 



^+f;M e rU k 



n=2 



vilH hj„=fc 7=1 



u BS (a ). 



Now, we insert expansions ([2]) and (J7J into ([5]) and collect terms of like order in e 
0(1) : u - u B5 (a ), 

0(e fe ) : u k 



oo ^ / n \ 

= a k d„u BS (a ) + ]T ^ E IK ^Vo), 



(6) 



(7) 



fc > 1. 
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Solving the above equations for {ofc}^ we find 
0(f) : Co = ao, 

«K«*)= <"■= ». ,U 1 >:^l >: II". I.';;.. 1 "'"..)]. i i. .». 



d a u BS 

Remark 10. The right hand side of ((HJ involves only <jj for j < k — 1. Thus, the {cfcj-J^ can be found 
recursively. 

Remark 11. Note that d™u BS (a) is easily computed using 

5> ss (a) = JdX (^e*tf s <">) (i/'a j h)ipx ■ 

Explicitly up to 0(e 4 ) we have 

0(e) : o"i Ul 



, 2 x M 2 - 2Tgld>0 

U(e ) : cr 2 = , 

OctUq 

3^ W 3 - (CT 2 CTl<92 + ^cr^Mo 

U(£ J : cr 3 = ■ . 

OaU 

0(e 4 ) : ct 4 = 



_ Ui - {a z a x dl + \a\dl + \a 2 a\dl + ^cr 4 9 4 )u 



We summarize our implied volatility result in the following theorem: 

Theorem 12 (Implied Volatility). The implied volatility a £ defined in (0 is given explicitly by ^ where 
Co = a and {<Jk\'k' = i are given by ([5]). 

Remark 13. Everything we have done so far is exact . The accuracy of the implied volatility expansion © 
is limited only by the number of terms one wishes to compute. 

5 Example: CEV 



The constant elasticity of variance (CEV) model of 



Coxl (|I 9751 ) improves upon the Black-Scholes model by 



allowing the volatility to depend on the present level of the underlying through a local volatility function (in 
the variable y = logx) of the form a(y) = y/e aie /3y / 2 . This model has enjoyed wide success because (i) it 
allows for closed- form solutions for European option prices and (ii) when (3 < 0, the local volatility function 
increases as y — > — oo (x — > 0), which is consistent with the leverage effect and results in a negative implied 
volatility skew. Still, the CEV model has some shortcomings. First, volatility drops to zero as y — > oo 
(x — > oo). Second, the model does not allow the underlying to jump. 
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We can correct both of these shortcomings in our framework by allowing rj(y) = ep(y) := e@ v . In this 
setting, the dynamics of Y become 

dY t =a(Y t )dt+ (al + ea\e f)Yt \ 1/2 dW t + f z dN t (Y t - , dz) , Y = y, 

a(Y t ) - (co + e cie w ) - \ [a\ + e a\e^) - J v{Y t _ , dz) (e z -l-z), 

C = inf jt > : (co + s ae pY ') ds = £ j , 

v(Y t -,dz) = v (dz) +eeP Yt -vi(dz) 

We call this model the "CEV" or "Constant Levylasticity of Variance" model, as it combines aspects of both 
the CEV model and exponential Levy models. We point out some of the main features of this model below: 

1 /2 

• For /3 < 0, the local volatility function a(y) — (<Zg + eafe^ v ) behaves asymptotically like the CEV 
model a{y) ~ y/eaie^ v ' 2 as y — >• — oo (a; — > 0), reflecting the fact that volatility tends to increase 
as the asset price X drops (the leverage effect). As y — > oo (x — ► oo) the volatility function behaves 
like a constant ~ do- This is in contrast to the CEV model, in which volatility function goes 
unrealistically close to zero as y — > oo (x — > oo). 

• Jumps of size dz arrive as a Poisson process with an instantaneous state-dependent intensity of 
v(Yt-,dz). The state-dependent Levy measure v(y : dz) = vo(dz) + ee^ v v\{dz) behaves asymptoti- 
cally as u(y,dz) ~ ee^v\(dz) as y — > — oo and asymptotically like v{y 1 dz) ~ ^o(rf^) as j/ — > oo. Thus, 
both the jump intensity and jump distribution can change drastically depending on the value of y. 

• The process Y is killed with a state-dependent intensity of k(Y t ) = (co + £C\e^ Yt ). For (3 < the 
function fc(y) behaves asymptotically like ec\e^ v as y — > — oo, reflecting the fact that a default is more 
likely to occur as Y drops. As y — > oo the function k(y) be haves asymptotically as a constant k{y) <~ cq. 



This particular form of killing rate was first suggested by 



Carr and Linetskvl (2006). 



• Although both the volatility o~(y) and jump intensity v(y, dz) are unbounded as y — ¥ — oo (for j3 < 0), 
negative infinity is still a natural boundary for the process Y. This is due to the fact that as y — > — oo 
the killing rate k(y) is also unbounded. 

Remark 14. Note that ep ^ Cq°(R). Thus, strictly speaking, the CEV model does not belong to the class 
of models described in section[2j Nevertheless, for any yo > — oo the function ep is an element of C^°(Mo), 
where Mo := (yo,°o)- Thus, we expect the results of this section to remain valid for all y > yo for e yo small 
enough. The numerical experiments provided at the end of the paper support this notion. 
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We wish to find an expression for u n when r] = Bp. Noting that 

(ip^epipx) =5(X- fi- i0), 
we see that the (n + l)-fold integral Q collapses into a single integral 



in= I d\ V ™ — TT X\-ikf) (ipX, h) tpX-i n 

1 1 ^ lljy fc (<PA-ife/3 -(p\-ijp) I 1 Xi 1 



f'n-l 



e„/3 / dX V „„ , r TT XA-ifc/3 {ip\,h)il)x- (9) 



Remark 15. Although we have written the option price as an infinite series ©, from a practical standpoint, 
one may only compute u £ f=a u^ N ^ :— Xm=o £ n u n . For any finite N we may pass the sum through the integral 
appearing in @. Thus, for the purposes of computation, the best way express the approximate option price 
is 



dx(^h)^Y, £ne ^[T,Tv^ — 2 — T n 



XA-ifc/3 



Note, to obtain the approximate value of u £ , only a single integration is required. This makes our pricing 
formula as efficient as other models in which option prices are expressed as a Fourier-type integral (e.g. Levy 
processes, Heston model, etc.). 

Numerical Results 

Define the transition density p £ (t, y; yo) and the 0(s n ) approximation of the transition density p^ n \t, y; j/o) 

n 

p £ (t 7 y;y ) = E yo S y (Y t ), p {n \t,y;y ) = ^e k p k (t,y;y [) ). 

k=0 

In figures Q] and [2] we plot the approximate transition density p( n ' for a specific model within the CEV 
framework. In figure [3] we plot the implied volatility smile induced by a specific model within the CFV 
framework (see equation ([5]) for a definition of implied volatility). The smile is generated by computing u £ 
first (using Theorem [5]) and then by inverting u B (a £ ) numerically to obtain a £ . 

Next, for a particular model in which v§ = and cq — (the assumption of section [4j define the 0(e") 
approximation of the implied volatility 



k=0 



In figure 3] we provide a numerical example illustrating convergence of c' n ) to a £ . We compute a*-™- 1 using 
Theorem [12] Implied volatility is plotted as a function of the log-moneyness to maturity ratio, LMMR := 
{k — y)/t. Convergence is fastest for values of k near y and slows as k moves away from y. 
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6 Conclusion 



In this paper we introduce a class of local Ley models. Within our modeling framework, we obtain a formula 
(written as an infinite series) for the price of any European option. Additionally, we obtain an explicit 
expression for the implied volatility smile induced by a certain sub-class of our models. As an example of 
our framework, we introduce a CEV-like local Levy model, which corrects some short-comings of the CEV 
model; namely (i) our choice of volatility function does not drop to zero as the value of the underlying 
increases and (ii) our model permits the underlying asset to experience jumps. In the CEV-like framework, 
we show that option prices can be computed with the same level of efficiency as other models in which option 
prices are computed as Fourier-type integrals. 

Thanks 

The author would like to thank Bjorn Birnir for his helpful comments. 
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A Spectral theory of normal operators in a Hilbert space 



In this appendix we summarize the theory of normal o perators acting; on a H ilbert 



Reed and Simonl (|1980f ) and 



space. A d etailed 



Rudinl ( 1973 ) 



exposition on this topic (including proofs) can be found in 

Let !K be a Hilbert space with inner product (•, •). A linear operator is a pair (dom(A)^A) where dom(A) 
is a linear subset of J£ and A is a linear map A : dom(A) — > 5£. The adjoint of an operator A is an operator 
A* such that (Af,g) = (f,A*g),Vf E dom(A),g € dom(A*), where 

dom(A*) :={ 9 €J(:3fteJ{ such that (Af, g) = (/, h) V/ € dom(yi)}. 

An operator (dom(./l), .A) is said to be self-adjoint in 3-C if 

dom(A) = dom(.A* ) , (A/, g) = (f,Ag) V /, <? e dom(A) . 

Throughout this appendix, for any self-adjoint opera tor A, we will assume that dom(A) is a dense subset of JC. 



A densely defined self-adjoint operator is closed (see 



Rudinl (|1973f ). Theorem 13.9). An operator (dom(A),.A) 



is said to be normal in Jt if it is closed, densely defined and commutes with its adjoint: A* A = AA* . Clearly, 
every self-adjoint operator is a normal operator. 

Given a linear operator A, the resolvent set p(A) is defined as the set of A € C such that the mapping 
(A— Id A) is one-to-one and R\ := (A — Id A) -1 is continuous with dom(i?>) = !K. The operator R\ : — >• 5£ 
is called the resolvent. The spectrum a (A) of an operator A is defined as a(A) := C \ p{A). We say that 
A € cr(A) is an eigenvalue of .A if there exists iji G dom(A) such that the eigenvalue equation is satisfied 



(10) 



A function -0 that solves (|10p is called an eigenfunction of .A corresponding to A. The multiplicity of an 
eigenvalue A is the number of linearly independent eigenfunctions for which equation (|10[) is satisfied. The 
spectrum of an operator A can be decomposed into two disjoint sets called the discrete and essential 
spectra: cr(A) = ad(A) U a e (A). For a normal operator A, a number A € C belongs too^L/l) if and only if A 



Rudinl (|l973l) . Theorem 12.29). 



is an isolated point of a(A) and A is an eigenvalue of finite multiplicity (see 

A projection-valued measure on the measure space (C, 23(C)) is a family of bounded linear operators 
{E(B),B £ 3(C)} in J£ that satisfies: 

1. E(0) = and E(C) = Id. 

2. is an orthogonal projection. That is, E 2 (B) = E(B) and E(B) is self-adjoint: E*(B) = E(B). 



1 The essential spectrum may be further decomposed into the continuous spec trum and the residual spectrum. It can be 
shown that the residual spectrum of an ordinary differential operator is empty (see 



Roach 



(1982), page 184). 
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3. E(Af]B) = E(A)E{B). 

4. If B = B, and B, n ■ = for i ^ j then E{B) = lim JWOC £™ =1 E ( B j)> where the limit is in the 
strong operator topology. 

5. For every in /, g £ "K the set function fj,t !g (B) := (f,E(B)g) is a complex measure on 23(C). 

Theorem 16 (Spectral Representation Theorem). There is a one-to-one correspondence between normal 
operators A and projection-valued measures E on "K, the correspondence being given by 



A= / XE(dX). 



If * s a B orel function on C then 



g(A) = I g(X)E(dX), 



Proof. See 



domfoGA)) = {/ € JC : f \g(X)\ 2 Nj (dX) < oc}. (11) 

J<j(A) 



Rudinl (|l973l ) Theorems 12.21 and 13.33. 



□ 



As a practical matter, if A is an differential operator acting on a Hilbcrt space L 2 (M.,dy), then the operators 
defined by (|ll|l can be constructed by solving the proper and improper [f| eigenvalue problems 



proper: 
improper: 



Aip n = <j> n ijj n , 
A ip\ — cj)\ tpx, 



4>n € CTd(~4), 

4>x e a e {A) 1 



For the improper eigenvalue problem one extends the domain of A to include functions all functions tp for 
which Af makes sense and for which the following boundedness conditions are satisfied 

lim \%p{y)\ 2 < oo. 

After normalizing, the proper and improper eigenfunctions A satisfy the following orthogonality relations 



The operator g(A) in (fTTj) is constructed as follows (see lHanson and Yakovlevl (|2002h . section 5.3.2) 



g(A)f= V g(<f>x)(ipx,f)i>x+ [ g{<t>x){^xJ)i>\dX 



\€tr d (A) 



It is not always easy to evaluate divergent integrals of the form (ipXi^X 1 ) an d verify that they are in fact 
delta functions S(X — X') . A method for d irectly obtaining properly normalised improper eigenfunctions can 
be found on page 238 of iFriedmanl ()1956l ). 



2 The term "improper" is used because the improper eigenvalues A ^ <?d{A) and the improper eigenfunctions ipx t since 
(i>X,">P\) = oo. 
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B Proof of Theorem [2] 



Our strategy is to show that A £ = Ao + enAi generates a semigroup CPf = exp(tA £ ). This will guarantee that 
u e (t,y) = y e t h(y) is an analytic function of e, which in turn, justifies the use of expansion Throughout 
this section we will work on the Hilbert space Di = L 2 (R, dy). We let do m(>/lj) = Cn°(R) and we note that 



(R) is a dense subset of Jf. Our analysis begins with a Theorem from 



Chernoffl ( 1972 ) 



Theorem 17. Let A be the generator of a Co contraction semigroup CP" = exp(L4) on a Banach space. Let 
e 23 be a dissipative operator with a densely defined adjoint. Assume that the inequality 

||eSu|| < a||u|| +b\\Au\\ , Vu € dom(A) 

holds for some a > and b < 1 (i.e., the operator sH is A-bounded with bound b < 1). Then the closure of 
A e := A + eH generates a C'o contraction semigroup CPf = exp(i.A e ). 

Remark 18. The operator Ao is the generator of a Cq contraction semigroup CP" = exp(L4o) on !K. 
To show that er\A\ is dissipative, the following Theorem will be useful: 

Theorem 19. Let A be a linear operator with domain dom(yi) = Cq°(M). Then A satisfies the positive 
maximum principle if and only if 

1 O , , „o . , , „ f , . , / g 



A = 2« 2 (y)9 2 + Ky)d + J v{y, dz) {e zd - 1 - l {z<R] zd) - c(y), (12) 
for some a{x) > 0, b(x) € K, c(x) >0,R€ [0, oo] and v{y,dz) satisfying 

v{y,dz) (l A z 2 ) < oo. 
Operators of the form (|12j) are called Levy- type operators. 

Proof. See Theorem 2.12 of|HobJ |l998l ). □ 



Rem ark 20. An operator that satisfies the positive maximum principle is dissipative (see 



Ethier and Kurtz 



(|1986l ). Lemma 4.2.1 on page 165). 



Remark 21. The operator er/Ai is clearly of the form (IT21 . Therefore, er\A\ is dissipative. 

Remark 22. The adjoint of er/Ai, given by (eijAi)* = eA\r), has domain doTtv{eA\r]) = Cq°(M), and is 
therefore densely defined in "K. 

Since .Ao generates a Co semigroup and er\A\ is dissipative and has a densely defined adjoint, we have only 
to show that et]A\ is .Ao-bounded with b < 1. 
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Proposition 23. Suppose rj G C^°(K) and e satisfies 

e < ini = , 

/or some a > and b < 1 (which is the condition given in Theorem^). Then er\A\ is Aa-bounded. 
Proof. First, observe that for any u £ C^(M.) we have 

\\A Q u\\ 2 = / dX\(^,u)\ 2 \^ x \ 2 , \\A lU \\ 2 = [ dA|(^, W )| 2 |xA| 2 , |M| 2 = / d\\[^u)f 



Thus, 



" AeR NI 2 -Ixa| 2 " NI 2 -Ixa| 2 

«■ O<a 2 + 6 2 |0 A | 2 - £ 2 ||r7|| 2 .|xA| 2 VA 

=* 0< |dA|(^A^)| 2 ( a 2 + 6 2 |0A| 2 -e 2 hl| 2 -|XA|- 

= a 2 || U || 2 + 6 2 ||^|| 2 -e 2 ||r/|| 2 -||yiiM|| 2 
<S> e 2 ||?/j| 2 • Hiull 2 < a 2 \\u\\ 2 + b 2 \\A u\\ 2 
||£??^iu|| 2 < a 2 \\u\\ 2 + b 2 \\A u\\ 2 
=> \\eri Aiu\\ < a \\u\\ + b \\A u\\ . 

where we have used ||?7^liu|| < \\r)\\ ■ ||.Aiu|| in the second-to-last step. 
The proof of Theorem [2] is complete. 



15 



n = 

1.0 

A 

1 \ o.s 




n = 1 




n = 2 

1.0 

A 

/ \ « 




1 > Q.6 

I I 

/ \ 




/ 1 1 1 <*6 




i I 

| A 06 

A ; /I 

/I 1 l\ 

/ A / \\ 

y/\ \ 




-3 -2 -1 1 




-3 -2-10 1 


n = 3 

1.0 

A 
t » 

1 l o.s 
i 1 

'/ \ 1 0.6 

J \ - 
// \\ 




n = 4 

1.0 

A 
i ' 

I l , o.s 

I \ 
I ■ 

A 1 06 

/' \ \ 
/' \ \ 

/' I Hi 

J I \ \ Utt 




n = 5 

1.0 

A 
i \ 

1 l , o.s 
i 1 

/ \ 16 

fi\ \ 
1 1 \ A I 

/ ' ^ \ t 4 
/ ' \\ 

j i y 
' >, 




-3 -: -i (i i 


-3 -2 -1 1 


-3 -2-10 1 


n = 6 

1.0 

,'\ 

I \ 
f 1 

1 1 08 
1 1 




n = 7 

1 1 o.s 

i ' 
i j 
i 1 




n = 8 

1.0 

/ * 

i l 08 
i 1 




| I Q.6 

n \ • ,, j 




| 1 06 
1 1 

A \ i 
/ ' \ \ 




1 ' ( °-6 

/' \ i 0.4 
/ 1 \ 1 




-3 -2 -1 (I 1 


-3 -2 -1 I) 1 


-3 -2-10 1 


n{dz)= — exp 





Figure 1: A plot of the approximate transition density p^ n '(t, y; 0) (solid blue) for different values of n. For 
comparison, we also plot (dashed black). Note that the density of Y t has a fat tail to the left, which is 
expected since the volatility and jump intensity increase as y — > — oo. The following parameters are used in 
these plots: a = 0.20, at = 0.10, c = 0.0, a = 0.001, s = 1.0, s x = 0.58, e = 1, /3 = -0.35, t = 1.0. 



1G 




Figure 2: For different values of n we plot of the approximate transition densities p( n ) (t, y\ yo) for yo = —0.6 
(solid red), yo — 0.0) (dashed black) and yo = 0.6) (dot-dashed blue). Note that as yo moves closer to — oo 
(i) the transition density becomes wider (ii) convergence of the approximate density p( n ) to the exact density 
p e requires more terms (iii) the densities have fatter tails on the left than on the right. All three phenomena 
are due to the fact that volatility and jump-intensity rise as y — > — oo. The following parameters are used 
in these plots: a = 0.20, a x = 0.10, c = 0.0, a = 0.001, s = 0.15, m = -0.10, s x = 0.15, m x = -0.10, 
£ = 1, -0.95, t = 1.0. 
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Figure 3: A plot of the implied volatility smile induced by a model with Gaussian jumps. Units on the 
horizontal axis are LMMR. Note that the steep skew, which is due to the fact that volatility and jump 
intensity increase as the value of the underlying drops. The following parameters are used in this plot: 
a = 0.20, ai = 0.10, c = 0.0, ci = 0.001, s = 0.2, s x = 0.1, m = -0.2, mi = -0.1, e = 1, /3 = -1.25, 
t = 0.50. 
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Figure 4: We plot cr^ n ' (solid blue) and a e (dashed black) as a function of LMMR. The following parameters 
are used in these plots: a = 0.25, a\ — 0.10, c = 0.00, c\ — 0.001, s\ — 0.15, mi = —0.10, e = 1, 
P = -1.25, t = 0.55, y = -0.01. 
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